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True ground states are evaluated for a 2d Ising model with 
random near neighbor interactions and ferromagnetic second 
neighbor interactions (the Randomly Coupled Ferromagnet). 
The spin glass stiffness exponent is positive when the abso- 
lute value of the random interaction is weaker than the fer- 
romagnetic interaction. This result demonstrates that in this 
parameter domain the spin glass like ordering temperature 
is non-zero for these systems, in strong contrast to the 2d 
Ed wards- Anderson spin glass. 
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INTRODUCTION 

For more than two decades the intensive numerical 
work on the spin glass (SG) problem has been concen- 
trated almost exclusively on the Edwards- Anderson Ising 
spin glass : Ising spins on a regular [hyper] cubic lat- 
tice with random near neighbor Gaussian or binomial 
interaction distributions [jjj. The many possible alter- 
native Ising systems with a randomness ingredient have 
hardly been touched on and such results as exist have 
been largely ignored. 

One such family of alternative systems was proposed 
in g|. It consists of a 2d square lattice of L x L Ising 
spins <7j = ±1 with uniform ferromagnetic second near 
neighbor interactions of strength J, plus random near 
neighbor interactions Jij = ±AJ; we will refer to it as 
the RCF (Randomly Coupled Ferromagnet) model ||. 
It is described by the following Hamiltonian: 



Z2 J(T i a j 



(1) 



For each realization of the randomness the Jij are drawn 



with the constraint Jij = to reduce fluctuations. 

As the spins are coupled through the ferromagnetic 
second near neighbor interactions, the system can be 
partitioned into two inter-penetrating sublattices in 
checkerboard-like fashion. In the limit A = the two sub- 
lattices order ferromagnetically and independently below 
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the Onsager temperature T — 2.27 J. Because each sub- 
lattice can order up or down, there are four degenerate 
ground states. As was pointed out in for non-zero A 
the near neighbor interactions can be considered in terms 
of effective random fields exerted by each sublattice on 
the other, so that for finite A and large enough L the 
ferromagnetic sublattice ordering is expected to be bro- 
ken up, as in the 2d random field Ising (RFI) model [||. 
The ground state will consist of coexisting domains of 
each of the four types : up/up, up/down, down/up and 
down/down. The question is: is the break-up accompa- 
nied by paramagnetic order down to T = ? 

A number of Monte Carlo simulations were performed, 
and it was concluded on the basis of standard numerical 
criteria (2|^] that when the ratio A is less than about 1, 
the RCF systems show spin glass like ordering at a fi- 
nite temperature, whereas 2d Edwards Anderson ISGs 
are paramagnetic down to T — [0. The finite order- 
ing temperature interpretation was strongly questioned 
by Parisi et al js| who criticized the initial work || on 
the grounds that the results were restricted to relatively 
small sample sizes L. On the basis of Monte Carlo data 
obtained on rather larger samples Parisi et al suggested 
that the RCF systems are always paramagnetic down to 
T = 0, like the Edwards Anderson ISGs. Further large 
sample Monte carlo results ||] however indicated finite- 
temperature ordering. 

Here we present data from ground-state configuration 
evaluations which show unambiguously that RCF sys- 
tems indeed exhibit finite temperature SG like ordering 
for A less than about 1. This opens up new and intrigu- 
ing possibilities for the testing of fundamental properties 
of complex ordered systems at finite temperatures in a 
2d context. 



ALGORITHM 

In the present work, ground state configurations have 
been found for periodic boundary conditions, and for the 
case where in one direction the boundary conditions are 
switched to anti-periodic. By comparing the ground state 
energies of the different boundary conditions for each 
realization conclusions on the ordering behavior can be 
obtained. Similar studies were performed for simple d- 
dimensional EA spin glasses in d ~ 2 0, d = 3 0, and 
d = 4§. 

For readers not familiar with the calculation of spin- 
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glass ground states now a short introduction to the sub- 
ject and a description of the algorithm used here are 
given. A detailed overview can be found in 0] 

The concept of frustration |ic| l is important for under- 
standing the behavior of ± J Ising spin glasses. The sim- 
plest example of a frustrated system is a triple of spins 
where all pairs are connected by antiferromagnetic bonds, 
see fig. [I]. A bond is called satisfied if it contributes with a 
negative value to the total energy by choosing the values 
of its adjacent spins properly. For the triangle it is not 
possible to find a spin-configuration were all bonds are 
satisfied. In general a system is frustrated if closed loops 
of bonds exists, where the product of these bond- values 
is negative. For square and cubic systems the smallest 
closed loops consist of four bonds. They are called (ele- 
mentary) plaquettes. 

As we will see later the presence of frustration makes 
the calculation of exact ground states of such systems 
computationally hard. Only for the special case of the 
two-dimensional system with periodic boundary condi- 
tions in no more than one direction and without external 
field a polynomial-time algorithm is known ID]]. In gen- 
eral only methods with exponential running times are 
known, on says the problem is NP-hard Jl2]. Now for 
the general case three basic methods are briefly reviewed 
and the largest system sizes which can be treated are 
given for three-dimensional systems, the standard spin- 
glass model, were data for comparison is available. 

The simplest method works by enumerating all 2 N 
possible states and has obviously an exponential run- 
ning time. Even a system size of 4 3 is too large. The 
basic idea of the so called Branch- and- Bound algorithm 
p3[ is to exclude the parts of the state space, where no 
low-lying states can be found, so that the complete low- 
energy landscape of systems of size 4 3 can be calculated 

A more sophisticated method called Branch- and- Cut 
^5| , ^6| works by rewriting the quadratic energy function 
as a linear function with an additional set of inequali- 
ties which must hold for the feasible solutions. Since not 
all inequalities are known a priori the method iteratively 
solves the linear problem, looks for inequalities which are 
violated, and adds them to the set until the solution is 
found. Since the number of inequalities grows exponen- 
tially with the system size the same holds for the com- 
putation time of the algorithm. With Branch-and-Cut 
anyway small systems up to 8 3 are feasible. 

The method used here is able to calculate true ground 
states @ up to size 14 3 . For two-dimensional systems, 
as considered in this paper, sizes up to 50 2 can be 
treated. The method is based on a special genetic algo- 
rithm fl7| , |l8|| and on Cluster-Exact Approximation [ p^| . 
CEA is an optimization method designed specially for 
spin glasses. Its basic idea is to transform the spin glass 
in a way that graph-theoretical methods can be applied, 
which work only for systems exhibiting no frustrations. 
Next a description of the genetic CEA is given. 

Genetic algorithms are biologically motivated. An op- 



timal solution is found by treating many instances of the 
problem in parallel, keeping only better instances and 
replacing bad ones by new ones (survival of the fittest) . 
The genetic algorithm starts with an initial population 
of Mi randomly initialized spin configurations (= indi- 
viduals) , which are linearly arranged using an array. The 
last one is also neighbor of the first one. Then n a x Mj 
times two neighbors from the population are taken (called 
parents) and two new configurations called offspring are 
created. For that purpose the triadic crossover is used 
which turned out to be very efficient for spin glasses: a 
mask is used which is a third randomly chosen (usually 
distant) member of the population with a fraction of 0.1 
of its spins reversed. In a first step the offspring are 
created as copies of the parents. Then those spins are se- 
lected, where the orientations of the first parent and the 
mask agree 20 . The values of these spins are swapped 



between the two offspring. Then a mutation with a rate 
of p m is applied to each offspring, i.e. a fraction p m of 
the spins is reversed. 

Next for both offspring the energy is reduced by ap- 
plying CEA: The method constructs iteratively and ran- 
domly a non-frustrated cluster of spins. Spins adjacent 
to many unsatisfied bonds are more likely to be added to 
the cluster. During the construction of the cluster a local 
gauge-transformation of the spin variables is applied so 
that all interactions between cluster spins become ferro- 
magnetic. 

Fig. ^ shows an example of how the construction of 
the cluster works using a small spin-glass system. For 2d 
± J spin glasses each cluster contains typically 70 percent 
of all spins. The non-cluster spins act like local magnetic 
fields on the cluster spins, so the ground state of the clus- 
ter is not trivial. Since the cluster has only ferromagnetic 
interactions, an energetic minimum state for its spins can 
be calculated in polynomial time by using graph theoreti- 
cal methods [^l]g2) : an equivalent network is constructed 
|23|, the maximum flow ||,|| is calculated f\ and the 
spins of the cluster are set to their orientations leading 
to a minimum in energy. This minimization step is per- 
formed n m i n times for each offspring. 

Afterwards each offspring is compared with one of its 
parents. The pairs are chosen in the way that the sum 
of the phenotypic differences between them is minimal. 
The phenotypic difference is defined here as the number 
of spins where the two configurations differ. Each parent 
is replaced if its energy is not lower (i.e. not better) 
than the corresponding offspring. After this whole step is 
done n Q x Mi times, the population is halved: From each 



1 Implementation details: We used Tarjan's wave algorithm 
together with the heuristic speed-ups of Traif. In the con- 
struction of the level graph we allowed not only edges (v,w) 
with level(«;) = level (u)+l, but also all edges (v,t) where t is 
the sink. For this measure, we observed an additional speed- 
up of roughly factor 2 for the systems we calculated. 
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pair of neighbors the configuration which has the higher 
energy is eliminated. If more than 4 individuals remain 
the process is continued otherwise it is stopped and the 
best individual is taken as result of the calculation. 

The representation in fig. ^ summarizes the algorithm. 

The whole algorithm is performed np times and all 
configurations which exhibit the lowest energy are stored, 
resulting in n g statistical independent ground state con- 
figurations. The running time of the algorithm with suit- 
able parameters chosen (see Table I) grows exponentially 
with the system size. On a 80Mhz PowerPC processor 
a typical L = 40 instance takes 3 hours (15 hours for 
L = 56). 



RESULTS 

In this work ground states of the RCF are studied for 
system sizes up to L = 56 and values of A — 0.5, 0.7, 0.9, 
and 1.1. Usually 1000 different realizations were treated, 
each submitted to periodic (pbc) and antiperiodic (apbc) 
boundary conditions in one direction and always pbc in 
the other direction. The apbc are realized by inverting 
one line of bonds in the system with pbc. Because of 
the enormous computational effort, for the largest sys- 
tem sizes only realizations with A = 0.7 where consid- 
ered with large statistics (and about 100 realizations with 
L = 56, A = 0.9). 

The periodic ground states give a direct measurement 
of the T = break up length Lb at each value of A, which 
is defined as follows: For small enough L the ground 
states will always be such that there is a full ferromag- 
netic ordering within each sublattice. With increasing L, 
more and more samples will be found with ground states 
having at least one of the sublattices incompletely fer- 
romagnetic. The break up length Lb is defined p^j as 
the value of L above which more than half the samples 
do not have pure ferromagnetic order in each sublattice. 
For the binomial RFI model, Lb ~ 5.5exp(2/A 2 ) where 
A is the strength of the random field [E6[ . For the RCF 
the values of Lb are shown against A - in Figure 1. It 
was suggested in || that by analogy with the RFI results 
|p6| Lb(\) could be expected to vary as exp(l/(4A 2 )). In 
fact the data points for the true ground states lie on 
the line Lb ~ 3.2 exp(0.62/A 2 ). For the particular cases 
A = 0.5 and A = 0.7, Lb « 45 and « 10 respectively. 
With the wisdom of hindsight, it can be seen that the 
measurements done in [^MpJ for A = 0.5 were mainly in 
the regime L < Lb while for A = 0.7 the larger samples 
were well in the regime L > Lb- 

A "typical" ground state for A = 0.7 and L — 56 is 
shown in Figure [| All four possible types of domains oc- 
cur. Because of the discrete structure of the interaction 
usually the ground state is degenerate. But in contrast 
to the EA spin glasses with only ± J near neighbor inter- 
actions, where a complex ground-state landscape exists, 
the structure of the degeneracy is trivial for A < 1: the 



whole system may be flipped, sometimes it is possible 
to flip both sublattices independently, and usually some 
small clusters occur with can take two orientations. 

But for studying whether the model exhibits long 
range order or not, it is sufficient to concentrate on 
the ground-state energies Ep,Eap for periodic and an- 
tiperiodic boundary conditions. The energy differences 
A = Ep — Eap give information about whether a sys- 
tem exhibits some kind of stiffness against perturbations 
of the boundary, i.e. about the presence of order. A is 
called the stiffness energy. For samples with the same set 
of interactions the stiffness can be analyzed in terms of 
the size dependence of the average (A) and of the width 
W = v/f7 2 (A) of the distribution P(A). For A = 0.7 
the distribution is presented in Fig. pi The inset shows 
the behavior of the average stiffness energy as a func- 
tion of L for all four values of A. For system sizes larger 
the breakup length and A > 0.7 the stiffness energy de- 
creases, indicating that no ferromagnetic long range order 
is present in the system. For A = 0.5 the breakup length 
is very large, so the asymptotic behavior is hardly visible, 
but (A) seems to fall for L > 28. From direct evaluation 
of the magnetization (see Fig and Fig. ||) we conclude 
that no ferromagnetic order should be present beyond an 
upper limit A = 0.27(8). For smaller values of A noth- 
ing can be concluded from our data. Furthermore, for 
smaller values of A it remains possible that the ground 
states of the RCF model do not exhibit ferromagnetic 
ordering ordering for any finite value of the relative cou- 
pling constant A. 

From standard relationships |2^ , p8|] one can write 
(A) ~ L 6p with dp the ferromagnetic stiffness exponent, 
and W ~ L 6sG with Osg the spin glass stiffness exponent. 
Positive values of the exponents indicate a long range or- 
der. Because of the small system sizes an evaluation of 
the ferromagnetic exponent is difficult. From the results 
presented in Fig. ^ we find an asymptotic (L — > oo) value 
of F = -2 (A = 0.9,1.1). 

Now we turn to the question whether some kind of 
spin-glass order is present in the system. This can be 
investigated by analyzing The dependence of the variance 
cr 2 (A) of the stiffness-energy distributions on the system 
size, the result is shown in Fig. ^| For small system sizes 
the variance grows for all values of the coupling constant 
A. In order to exclude finite-size effects, only systems 
larger than the breakup length Lb{\) should be taken into 
account. Above Lb there is a good linear size dependence 
of \ogW{L) against logL, with SG = 0.59(8), 0.29(1), 
0.09(5), and -0.16(2) respectively for A = 0.5,0.7,0.9, 
and 1.1. 

The values of Osg against A are shown in the inset of 
Figure ^. The result for A = 0.5 is not very reliable, be- 
cause the largest system size is of the order of the breakup 
length. In the log-log plot the datapoints for A = 0.5 ex- 
hibit a negative curvature, thus the asymptotic value of 
Osg m ay be smaller than 0.59. For the other systems the 
breakup length is quite small, so the results give unam- 
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biguous evidence for spin glass like ordering in the large 
size limit, with a non-zero ordering temperature. Espe- 
cially for A = 0.7, where w 10, the result a (A) > 
is very reliable. Thus, it is indeed not necessary to carry 
out further calculations with larger systems to prove the 
fact, there there are values of the coupling constant giv- 
ing rise to an ordered spin glass phase in the RCF. The 
limiting value A c above which 6sg is negative is very close 
to 1.0; A c would correspond to the highest value at which 
the ordering temperature is non-zero, in good agreement 
with the initial estimate from the Monte Carlo work Q. 

CONCLUSION 

We have calculated ground states of the Randomly 
Coupled Ferromagnet for different values of the spin- 
glass coupling constant A and with periodic as well as 
antiperiodic boundary conditions. By using the genetic 
cluster-exact approximation algorithm, we were able to 
treat system sizes up to N = 56 x 56. The breakup length 
was calculated for each value of A. From the calculation 
of T = stiffness energy it could be concluded that be- 
low A c w 1 the RCF exhibits an ordered spin glass like 
phase at finite temperature. It should be stressed again 
that for A > 0.5 the largest system sizes are well beyond 
the breakup length, so no changes are to be expected for 
larger system sizes. For A < 0.5, especially if one likes to 
test whether the model exhibits ferromagnetic ordering, 
ground states calculations of larger systems are needed to 
study the behavior in more detail. Unfortunatley, these 
studies are beyond the power of current computers and 
algorithms. 

Although the zero-temperature stiffness exponent val- 
ues give no direct information on the ordering temper- 
atures, the present results are consistent with the the 
conclusions drawn in [^|J|] where Monte Carlo estimates 
of the critical temperatures were made using the finite- 
size scaling of the spin glass susceptibility and the form 
of the time dependence of the autocorrelation function 
relaxation. Ordering temperatures were estimated to be 
close to 2.0 for A = 0.5 and 0.7, dropping to zero near 
A = 1. Rather remarkably the T = crossover as a 
function of L at Lt, appears to have little effect on the 
behavior of the SG susceptibility as a function of size 
in the temperature region close to T g jl^]. However for 
A = 0.5 Parisi et al ||] observed weakly non-monotonic 
behavior of the Binder parameter with L for sizes that 
we now know to be in the region of the crossover. 

Since the existence of a spin glass like phase at finite 
temperature now has been established definitely, it would 
be instructive to carry out further careful Monte Carlo 
measurements for sample sizes well in the regime L > Li, 
and over a range of A values. Is the physics of the 2c? RCF 
above, at, and below the ordering temperature strictly 
analogous to that of the standard Edwards Anderson ISC 
at dimensions where there is finite temperature ordering? 



To what extent could the RCF enlighten us concerning 
problems which in the Edwards Anderson ISC context 
have remained conflictual for more than twenty years ? 
The fact that the RCF lives on a 2d lattice rather than 
in a higher dimension should facilitate understanding of 
the fundamental physics of ordering in complex systems. 

Finally, there may even be possible experimental re- 
alizations of systems where quasi-two-dimensional mag- 
nets form short range clusters with local ferromagnetic 
or antiferromagnetic order, with random frustrated in- 
teractions linking these clusters together. Examples of 
promising behaviour of this sort are Fe compound with 
halogens |^9| where it might be interesting to look at the 
data again in view of the present results. 
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TABLE I. Simulation parameters: L — system size, Mi = 
initial size of population, v = average number of offsprings per 
configuration, n m i n = number of CEA minimization steps per 
offspring, p m = mutation rate, n.R=number of independent 
runs per realization. 
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FIG. 1. The simplest frustrated system: a triple of spins, 
each pair of spins connected by antiferromagnetic bonds 
(dashed lines). It is not possible to satisfy all bonds. 




FIG. 2. Example of the Cluster-Exact Approximation 
method. A part of a spin glass is shown. The circles repre- 
sent lattice sites/spins. Straight lines represent ferromagnetic 
bonds the jagged lines antiferromagnetic interactions. The 
top part shows the initial situation. The construction starts 
with the spin at the center. The bottom part displays the fi- 
nal stage. The spins which belong to the cluster carry a plus 
or minus sign which indicates how each spin is transformed, 
so that only ferromagnetic interactions remain inside the clus- 
ter. All other spins cannot be added to the cluster because it 
is not possible to multiply them by ±1 to make all adjacent 
bonds positive. Please note that many other combinations of 
spins can be used to build a cluster without frustration. 
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algorithm genetic CEA({Jy}, M i; n OI p m , n min ) 
begin 

create Mi configurations randomly 

while (Mi > 4) do 

begin 

for i = 1 to n x Mi do 

begin 

select two neighbors 

create two offspring using triadic crossover 
do mutations with rate p m 
for both offspring do 
begin 

for j = 1 to n m j n do 

begin 

construct unfrustrated cluster of spins 
construct equivalent network 
calculate maximum flow 
construct minimum cut 
set new orientations of cluster spins 
end 

if offspring is not worse than related parent 
then 

replace parent with offspring 

end 
end 

half population; Mi = Mi/ 2 
end 

return one configuration with lowest energy 
end 

FIG. 3. Genetic Cluster-exact Approximation. 
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FIG. 5. Typical ground state of one RCF realization 
(L = 56) for A = 0.70 with periodic boundary conditions in all 
directions. Two different symbols (white/black square), (un- 
filled/filled diamond) are used to represent the orientations 
on the different sublattices. 
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FIG. 6. Distribution of stiffness energies A at A = 0.7 for 
different system sizes L = 14, 28, 56. The inset shows the 
mean stiffness energy (A) as a function of L for A = 0.5, 0.7, 
0.9, and 1.1. 



FIG. 4. Breakup length Lb as a function of 1/A 2 . The solid 
line shows a fit to the function Lb(X) = aexp( — b/X 2 ). 
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FIG. 7. Average magnetization m as a function of the 
strength A of the spin-glass couplings for different systems 
sizes 5 < L < 56. 
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FIG. 9. Width cr^(A) of the distribution of stiffness 
energies as a function of the system size L > 10 for 
A = 0.5, 0.7, 0.9, 1.1. The solid lines are fits to algebraic func- 
tions of the form cr 2 (L) = gL BsG . The inset shows the values 
of the exponent for different values of A. 




FIG. 8. Threshold value A c as a function of the system 
size L. The threshold value is determined from m(A c ) = 0.9. 
The solid line shows a fit to a function A C (L) = Aoo + eL~, 
resulting in A^ = 0.27(8), / = -0.53(10). 
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